arXiv: 1506.00587v2 [cond-mat.quant-gas] 26 Nov 2015 


Coherent driving and freezing of bosonic matter wave in an optical Lieb lattice 

Shintaro TaieQ Hideki Ozawa, Tomohiro Ichinose, Takuei Nishio, Shuta Nakajima, and Yoshiro Takahashi 
Department of Physics, Graduate School of Science, Kyoto University, Japan 606-8502 

(Dated: November 30, 2015) 

While kinetic energy of a massive particle generally has quadratic dependence on its momentum, 
a flat, dispersionless energy band is realized in crystals with specific lattice structures. Such macro¬ 
scopic degeneracy causes the emergence of localized eigenstates and has been a key concept in the 
context of itinerant ferromagnetism. Here we report the realization of a “Lieb lattice” configuration 
with an optical lattice, which has a flat energy band as the first excited state. Our optical lattice po¬ 
tential possesses various degrees of freedom about its manipulation, which enables coherent transfer 
of a Bose-Einstein condensate into the flat band. In addition to measuring lifetime of the flat band 
population for different tight-binding parameters, we investigate the inter-sublattice dynamics of the 
system by projecting the sublattice population onto the band population. This measurement clearly 
shows the formation of the localized state with the specific sublattice decoupled in the flat band, 
and even detects the presence of flat-band breaking perturbations, resulting in the delocalization. 

Our results will open up the possibilities of exploring physics of flat band with a highly controllable 
quantum system. 

PACS numbers: 34.50.-s, 67.85.-d 


Many-body properties of a quantum system show dras¬ 
tic changes according to the geometry of an underlying 
lattice structure. One of the textbook examples is an an- 
tiferromagnet on a frustrated lattice geometry [lj> where 
the geometric frustration prevents spins from Neel order¬ 
ing and the system exhibits more nontrivial, correlated 
ground state. Specific lattice geometry can also induce 
frustration of kinetic energy. In this case a macroscopic 
number of momentum eigenstates are energetically de¬ 
generate, forming a dispersionless flat band. Flat bands 
have been playing an important role in theoretical study 
of itinerant ferromagnetism, as the presence of interac¬ 
tion lifts the bulk deg eneracy and chooses the ferromag¬ 
netic ground state [2j-|4|. 

Ultracold atomic gases in a periodic potential gener¬ 
ated by standing waves of laser light (optical lattices) 
have made a great success in realizing controllable many- 
body systems described by well-defined theoretical mod¬ 
els of interest, such as the Hubbard model H, @. Re¬ 
cently, increasing experimental efforts have been made 
to create and investigate non-standard (other than sim¬ 
ple cubic) optical lattices which have unique geometric 
features [3-E3- Among them, a kagome lattice 0 is 
a promising candidate for exploring flat band physics. 
However, the negative hopping amplitude, which is nat¬ 
ural for ultracold atoms, forces the flat band to be ener¬ 
getically highest among the three s-orbitals, and so far 
phenomena characteristic to the flat band have not been 
reported. A different type of lattice structure known as a 
Lieb lattice, also referred as a decorated square lattice, or 
a d-p model in the analogy to CuCU plane of a High-T c su¬ 
perconductor [l(|, similarly has a flat band as the second 
(first excited) band. It consists of two sublattices: one of 
them forms a standard square lattice (the A-sublattice 
in Fig. QJi) and the other lies on every side of the square. 
For convenience, we further divide the latter into the B¬ 


and C- sublattices. The single particle energy spectrum 
in the tight binding limit (Fig. [ljr) has the characteristic 
flat band and the Dirac cone on the corner of the Brillouin 
zone. This lattice satisfies the criteria for the occurrence 
of Lieb’s ferrimagnetism, which states that the half-filled 
spin-1/2 fermions exhibits nonzero magnetization for a 
positive on-site interaction Q. Also for bosonic systems, 
a flat band gives a fascinating question of whether con¬ 
densation is possible in the presence of kinetic energy 
frustration. Theoretical investigation predicts supersolid 
order for a flat band 17]. The Lieb lattice was recently 
realized in a photonic lattice jTs]] where the electric field 
passing through the lattice emulates time evolution of 
single-particle wave function. Also, Bose condensation of 
polaritons in one-dimensional analogue of the Lieb lattice 
(sawtooth lattice) has also been reported |l9j] . Optical- 
lattice realization of the Lieb lattice provides powerful 
quantum simulation of a closed, interacting many-body 
system. 


In this paper, we present manipulation and detection of 
bosonic matter-wave in an optical Lieb lattice. An atomic 
condensate is coherently loaded into the flat band and its 
time evolution is measured in two ways: quasimonentum- 
and sublattice-resolved detection, which reveal the local¬ 
ized character of the flat band wave function. Relatively 
short lifetime of atoms in the 2nd band was observed, 
while it can be made longer by increasing the band gap 
to the lowest band. This work paves the way to ex¬ 
perimental study of flat band physics with cold atoms. 
Using Fermi gases with the Fermi energy lying at the 
flat band can avoid the lifetime problem and will pro¬ 
vide ideal playground for investigating flat band ferro¬ 
magnetism (20-1221] and topological phases with artificial 
gauge fields |23|. 
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FIG. 1. Optical Lieb lattice, a, Lieb lattice. A unit cell is 
indicated by the green square, b. Tight-binding energy band 
structure of the Lieb lattice, c, Experimental realization of 
the Lieb lattice. Black arrows indicate polarizations of the 
lattice beams, d, Lattice potential for (siong, Sshort, s diag) = 
(8,8, 9.5) at (j) x = 4>z = 0 and ip = 7t/2. e, Band structures 
of the optical Lieb lattice at: (si on g, Sshort, Sdiag) = (8, 8,9.5) 
(red dashed), (34,34,37.4) (black solid). 


FORMATION OF AN OPTICAL LIEB LATTICE 

We construct the optical Lieb lattice by superimposing 
three types of optical lattices (see Fig. [TJ:- and[]Jl), leading 
to the potential 

V ( x , z) = - cos 2 (fc L 2 ’) - 4ng cos 2 (k h z) 

~ Khort cos 2 (2fc L a: + <f> x ) - V^ It cos 2 (2 k L z + <p z ) 
- Fdiag cos 2 [fc L (x - z) +1p\ , (1) 

where z indicates the direction of gravity. Here, k l = 
2ir/\ is a wavenumber of a long lattice (with a depth 
Mong), for which we choose A = 1064 nrn. A short lattice 
(bshort) is formed by laser beams at 532 nrn. A diag¬ 
onal lattice (Vdi a g) with the wave number \/2fcL is re¬ 
alized by interference of the mutually orthogonal laser 


beams at 532 nm along the x- and ^-directions. Com¬ 
pared to the proposals nm, equation m lacks a di¬ 
agonal lattice with x + z spatial dependence. While this 
causes slightly larger discrepancy from the ideal tight- 
binding energy bands, the sufficiently deep lattices can 
reproduce the desired energy spectrum including the flat 
band, as seen in Fig. |TJ?. In the following, we specify 
each lattice depth in unit of long-lattice recoil energy as 

(siong; Sghort? Sdiag) = (Cong- Kdiort: Clia.g)/^R ■ where 

4 1064 > = 1i 2 k1/ (2m) and m is the atomic mass of 174 Yb. 
The motion of free bosons in the x-z plane is governed 
by tight-binding Hamiltonian 

F =(ajoj+H.c.) + Es^fn, (2) 

(id) S=A,B,C ieS 


where hi is the annihilation operator on site i and 
hi = ajdj. The nearest-neighbor hopping amplitude 
J is mainly determined by V'short, whereas the other 
lattice depths set an energy offset Es of each sublat¬ 
tice. Excepting the contribution from the zero-point 
energies of each potential well, they are approximately 
given by E A ~ -V^ - V^ g , E B ~ -V^ g ~ F diag , 

and Ec ~ — ^[ong — ^diag? which are able to be indepen¬ 
dently controlled by tuning the lattice depths. In the y- 
direction which is perpendicular to the Lieb lattice plane, 
the atoms are weakly confined in a harmonic trap (one¬ 
dimensional tube configuration), unless otherwise speci¬ 
fied. In addition to each lattice depth, three phase pa¬ 
rameters should be set to 4> x = <p z = 0 and 0 = 7r/2 in 
order to realize the characteristic three-sublattice struc¬ 
ture of the Lieb lattice depicted in Fig. [Hi (Methods). 


LOADING BEC INTO A FLAT BAND BY PHASE 
IMPRINTING 

One of the fundamental properties of flat bands is 
the localization of the wave function as a consequence 
of quantum-mechanical interference of traveling matter 
waves. The localization is due to a purely geometric ef¬ 
fect, as we briefly explain below. The Hilbert space for 
the Lieb lattice in the tight-binding regime is spanned by 
the quasimomentum eigenstates of each sublattice | k,S) 
(S = A,B,C). Nearest neighbor tunneling induces the 
momentum dependent coupling Jab = J cosk x d/2 be¬ 
tween the A- and B-sublattices, and similarly Jag = 
J cosk z d/2 between the A- and C-sublattices, where 
d = 532 nm is the lattice periodicity. The flat band states 
are the zero-energy eigenstates cos6*|fc,.B) — sin^lfejC) 
with tanf? = Jab/Jac, which have no amplitude on 
the A-sublattice. Consequently, a wave packet composed 
of the flat band states remains localized, as the tunnel¬ 
ing from a B-site and a C-site destructively interfere on 
the adjacent A-site. We explore this nature in the fol¬ 
lowing experiments. Here we note that the flat band in 
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FIG. 2. Coherent band transfer, a, Principle of the trans¬ 
ferring method, b, Absorption images reveal the coherent os¬ 
cillations between the |B) + |C) and \B) — \C) states. In the up¬ 
per left image, the first three Brollouin zones are displayed by 
white, green, and red lines, respectively, c, Oscillating behav¬ 
ior of the band population during phase imprinting in absence 
of lattice confinement along the {/-direction (blue circles), and 
with lattice confinement — V y cos 2 (2ki,y) (red squares). Solid 
lines are the fit results using the single particle solution of 
the Schrodinger equations (see main text). Error bars denote 
standard deviation of three independent measurements. 


the Lieb lattice is mathematically equivalent to “dark 
states” of laser-coupled A-type three level systems in 
atomic physics. Here, three sublattices correspond to the 
basis of three levels, tunneling amplitudes serve as laser- 
induced coupling, and energy difference of each sublattice 
plays a role of the detuning of laser. 

In a Lieb lattice, a flat band is realized as the first 
excited band, hence a BEC loaded adiabatically into an 
optical Lieb lattice is not populated in the flat band. 
However, tunability of our optical Lieb lattice enables 
to coherently transfer the population in the lowest band 
into the flat band by phase imprinting (see Fig. [2b)- The 
scheme is easily understood by considering tight-binding 


wave functions in each band. At zero quasimomentum 
and in the equal-offset condition Ea = Eb = Ec , a 
simple calculation gives 11st) = \A) + (\B) + \C))/\/2, 
12nd) = | B) - | C), and |3rd) = | A) - (|H) + \C))/y/2 
from the 1st to the 3rd band, where we omit the mo¬ 
mentum indices from the sublattice eigenstates (see also 
Supplementary Information). Taking advantage of rich 
controllability in our lattice potential, we can smoothly 
modify these eigenstates. With sufficiently large I4ia g 
(equivalently with large Ea — Eb,c ), the lowest Bloch 
state has essentially no amplitude in the A-sublattice, 
allowing to realize | B) + \C) state. Next, we apply sud- 
den change in one of the long lattice, This creates 

the energy difference between the B- and C-sublattices 
and the relative phase of the condensate wave function 
starts to evolve with a period 2irh/(Ec — Eb)- In the 
basis of the initial band structure, this time evolution is 
a coherent oscillation between 11st) and |2nd). 


The explicit procedure of the loading and detecting a 
condensate in the flat band is as follows. We adiabatically 
load a BEC of 2 x 10 4 174 Yb atoms into the Lieb lattice 
with (siong, Sshort, Sdiag) = (8, 8 , 20) and apply sudden in- 
crease of s lo „ to 26.4 for variable duration. At the same 
time, we ramp s s hort up to 20 to prevent tunneling dur¬ 
ing the band transfer. After this sequence, we return the 
lattice depths to the initial values and perform adiabatic 
turning off of the lattice potential in order to map quasi¬ 
momentum to free-particle momentum (band mapping) 


26i . [27]. Figure [5b shows the absorption images taken 


after 14 ms of the ballistic expansion, which reveal the 
interband dynamics of a condensate. At zero quasimo¬ 
mentum, atoms in the 2nd and 3rd band are mapped 
to the same point of the Brillouin zone. In addition, 
the finite spread of the condensate makes it difficult to 
precisely distinguish the population in the 2nd Brillouin 
zone from other neighboring zones. Therefore, Instead of 
plotting the population in the 2nd Brillouin zone, here 
we count atoms in the 1st Brillouin zone, and show the 
fraction of atoms in the other higher bands in Fig. [2b- For 
momentum space analysis, see Section III of Supplemen¬ 
tary Information. We fit the data with the function in the 
form a exp(— t/T)Fth{t)+b(l— exp(—t/r)), where Fth(t) is 
the numerical solution of the single particle Schrodinger 
equation. In fitting the data, we adopt sj OI [ during the 
band transfer as a free parameter, and obtain the best 
fit with s| oi j = 25.5, close to the expected value of 26.4. 
Although the oscillations involve non-negligible contribu¬ 
tions from the higher bands, at the half period of the first 
cycle we expect that > 75% of atoms are transferred to 
the 2nd band. This transfer method is also applicable in 
the presence of the lattice confinement along the y-axis, 
though the decay time of the oscillation r = 86(7) p,s is 
much shorter than the case of a weak harmonic confine¬ 
ment (lD-tube), t = 260(10) /zs. 
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RELAXATION DYNAMICS OF A FLAT BAND 

We measure the lifetime of atoms in the 2nd band of 
the optical Lieb lattice. After transferring to the 2nd 
band, we change the depth Sdiag of the diagonal lattice 
to control the energy gap between the 1st and 2nd band. 
As well as the band gap |28|, lifetime of a quantum gas in 
the excited band is strongly affected by the density over¬ 
lap with the states in the lower bands [291 . As we increase 
Sdiag, the average gap between the 1st and 2nd band be¬ 
comes smaller and at the same time their density profiles 
become similar to each other. In the opposite limit of 
shallow Sdiag, the band gap increases and two bands have 
no density overlap, as the lowest band mostly consists of 
the A-sublattice. We take a variable hold time in the lat¬ 
tice, followed by band mapping to count the atom num¬ 
ber in the excited bands. Typical absorption images are 
shown in Fig. [3j\. The decay curves displayed in Fig. [3jo 
show expected behavior of increasing lifetime with de¬ 
creasing Sdiag- In addition, increasing the gap makes the 
dynamics more clearly separate into two process: decay 
of the condensate within the 2nd band (middle image of 
Fig-EK) and decay of atoms into the lowest band (bottom 
image). We find that the curve is well fitted by double ex¬ 
ponential with the form aq exp(—t/ri)+a 2 exp(— t/r 2 )+b. 
The fast component t\ shows only weak dependence on 
Sdiag, whereas the slow component 72 shows over 20-fold 
changes from the smallest to largest Sdiag- We also ex¬ 
tract the lifetime of the condensate in the 2nd band, r c , 
by counting atoms on the corner of the 2nd Brillouin 
zone (Figure |3jr), and find similar behavior with t\. This 
implies that the intial fast decay is related to the decay 
of the condensate, which involves the decay to the lower 
band with faster time constant compared with the non- 
condensed atoms. 

LOCALIZATION OF A WAVE FUNCTION IN A 
FLAT BAND 

As described above, the most intriguing property of 
a flat band is the localization of the wave function at 
certain sublattice sites. In the case of the Lieb lattice, 
the wave function of the flat band vanishes on the A- 
sublattice. Here we reveal this property by observing the 
tunneling dynamics of a Bose gas initially condensed at 
the |—) = \B) — \C) state, and compare it to the dynamics 
of the state with opposite relative phase, |+) = \B) + \C). 
In order to observe real-space dynamics of the system, 
we perform projection measurement of the occupation 
number in each sublattice, which we call sublattice map¬ 
ping. In this method, first we quickly change the lattice 

potential to (( S long> S long)> S short, Sdiag) — ((8, 14), 20, 0). 
In this configuration, all three sublattices are energet¬ 
ically well separated from one another and the lowest 
three bands consist of the A-,B- and C-sublattice, re¬ 


spectively. This maps sublattice occupations to band 
occupations, which then can be measured by band map¬ 
ping technique. Figure |4ji shows the demonstration of 
this method, in the cases where atoms occupy only one 
of the sublattices. Note that the populations in the B- 
and C-sublattices are mapped to the 2nd Brillouin zones 
for the one-dimensional lattice along the x and z-axis, re¬ 
spectively. This is because the turning off of the diagonal 
lattice decouples these two directions and the fundamen¬ 
tal bands are labeled by the combination of band indices 
of ID lattices. 

We prepare the initial state |+) by simply loading a 
BEC into the Lieb lattice with deep Vdiag- On the other 
hand, the |—) state is obtained by applying band trans¬ 
fer method to the |+) state. Dynamics of these initial 
states after the lattice depths are changed to satisfy the 
equal-offset condition Ea = Eb = Eq is measured by 
the sublattice mapping. As shown in Fig. |4 )d, we re¬ 
veal qualitatively different behaviors of these two states: 
the |—) state shows a significant suppression of the A- 
sublattice occupancy, indicating the freezing of the tun¬ 
neling dynamics to the A-sublattice from the |—) state 
with only a slow decay to the A-sublattice, whereas the 
|+) state exhibits coherent oscillations between the A- 
and (HC)-sublattices. This clearly features the geomet¬ 
ric structure of the Lieb lattice mentioned above. Double 
exponential fit to the data for the |—) initial state yields 
T\ = 0.36(4) ms and r 2 = 5.5(9) ms, indicating that the 
leakage to the A-sublattice is caused by the decay to the 
lowest band. 

In the Bloch basis, the state |+) is expressed as 
11 st) — 13rd) and its time evolution is driven by the 
band gap AE^i which equals 4\/2 J in the tight bind¬ 
ing limit. After a half period nTi/AEs^ the state evolves 
to |A) = 11st) +13rd), leading to coherent tunneling to the 
A-sublattice. Similarly, it is possible to arrange the initial 
lattice depths so that the lowest Bloch state has the max¬ 
imum overlap with a certain superposition of 11st) and 
12nd). Sudden potential change to the Lieb lattice drives 
oscillation between the B- and C-sublattices, whose fre¬ 
quency gives the band gap A£ 2 ,i- We fit these data 
with a damped sinusoidal oscillation and compare the ex¬ 
tracted frequency with the result of single particle band 
calculations (see Fig. |4jc) - Qualitative behavior is well 
reproduced, while quantitative discrepancies are found. 
This is caused by interactions, as we present a system¬ 
atic study of the density dependence of the oscillation 
frequency in Section IV of Supplementary Information. 

We further investigate the tunneling dynamics of the 
|—) initial state by adding the perturbations which de¬ 
stroy the flatness of the second band. The flatness 
is robust against the independent change of nearest- 
neighbor tunneling amplitudes J x , J z along the x- and 
z-directions, and energy offset Ea, just as a dark state 
in a A-type three level system persists regardless of laser 
intensities and detuning from the excited state. How- 
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FIG. 3. Lifetime of atoms in the flat band, a, Absorption images for the lifetime measurement of the 2nd band with 
three different hold times, taken after 14 ms time-of-flight. The diagonal lattice depth is Sdiag = 9.5. The first three Brillouin 
zones are indicated by the white dashed lines. In the top image, the areas used to evaluate the lifetime of a condensate (r c ) 
are also displayed with the red squares, b, Decay of the flat band at (si on g, Sshort) = (8,8) and variable Sdiag- Solid lines are 
the fit results with double exponential curves. Error bars denote the standard deviation of three independent measurements, 
c, Lifetime of the flat band, n ,2 are the fast and slow decay time obtained from the data shown in b, respectively. r c is the 
e^ 1 lifetime of condensates. Error bars represent fitting error. 


ever, if the energy difference between the B- and C- 
sublattices is introduced - the two-photon Raman off- 
resonant case the flat band is destroyed. Note that the 
finite Eb — Ec induces population in the A-sublattice 
even at k = 0. On the other hand, the direct diago¬ 
nal tunneling between the B- and C-sublattices, which is 
another flat-band-breaking term existing in our system, 
keeps a dark state at k = 0 provided J x = J z . We cre¬ 
ate the energy difference by introducing the imbalance of 
Asiong = s^ g — £long ■ Figure 0J1 shows the time depen¬ 
dence of the A-sublattice population for the |—) initial 
state. It can be clearly seen that the coherent tunneling 
dynamics starts to grow as the lattice parameters deviate 
from the flat-band condition Asi ong = 0. 


METHODS 

Preparation of 174 Yb BEC 

After collecting about 10' atoms with a magneto¬ 
optical trap with the intercombination transition, the 
atoms are transferred to a crossed optical trap. Then 
we perform an evaporative cooling, resulting in an al¬ 
most pure BEC with about 10 5 atoms with no discern- 
able thermal component. 

All of the optical lattice experiments presented in this 
paper are subject to additional weak confinement due to 
a crossed optical dipole trap operating at 532 nm. Gaus¬ 


sian shape of laser beams for the trap and lattices impose 
a harmonic confinement on atoms, whose frequencies are 
(u] x >,u}y',u] z )/2ir = (147, 37,105) Hz at the lattice depths 
of (^long: a s hort:-^diag) — (8, 8,9.5). Here, the x- and y- 
axes are tilted from the lattice axes (x and y) by 45° in 
the same plane. 


Construction of optical Lieb lattice 

The relative phases between the long and short lat¬ 
tice (<fi x , <fi z ) can be adjusted by changing the frequency 
difference between these lattice beams [30|. The proper 
frequencies that realize the Lieb lattice (t f> x = (f> z = 0) 
are determined by analyzing the momentum distribution 
of a 174 Yb BEC released from the lattice, as in the case 
of the parameter if> of the diagonal lattice (Supplemen¬ 
tary Information). The relative phase between the long 
and short lattices at the position of atoms depends on 
the optical path lengths from common retro-reflection 
mirrors, and in general two phases <p x and <f> z are not 
equal. We shift the frequency of long lattice laser by an 
acousto-optic modulator (AOM) inserted in the path for 
the 0 -axis in order to simultaneously realize (fi x = cj) z = 0. 
Optimal frequency difference is sensitive to the alignment 
of the lattice beams and day-by-day calibration of the 
phases is needed. Typical drift of the required RF fre¬ 
quency for the compensation AOM is within 5 MHz. 

To stabilize the phase we construct a Michelson in- 
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FIG. 4. Tunneling dynamics in the Lieb lattice, a, Demonstrating the measurement of sublattice occupancy. Here, sub- 
lattice mapping technique is applied to atoms loaded into (left) ((s£2 g > s ion g )> Sshort, Sdia g ) = ((8, 8), 8, 0), (center) ((2, 8), 8, 19), 
and (right) ((8, 2), 8,19), corresponding to atoms in A-, B-, and C-sites, respectively, b, Measured tunneling dynamics of |+) 
and |— ) initial states in the Lieb lattice with (si ong , Sshort, Sdia g ) = (8, 8, 9.5). Solid lines are the fits to the experimental data 
with damped sinusoidal oscillation (for |+)) and double exponentials (for |—)). Inset shows dynamics of the |—) state for longer 
hold times. Error bars denote standard deviation. Illustration of tunneling process for each initial state is also shown on the 
right-hand side, c, Frequencies of coherent inter-site oscillations. Solid lines are the calculated band gap between the 1st and 
2nd (red) and the 1st and 3rd bands (blue). Error bars denote fitting error, d, Bending flat band. Dynamics of the |—) state in 
the presence of imbalance Asi ong = s^ g — shows restoration of coherent dynamics. Error bars denote standard deviation. 


terferometer along the optical path of the diagonal lattice 
with frequency stabilized 507 nm laser. The interferom¬ 
eter has two piezo electric transducer (PZT)-mountded 
mirrors one of which is shared with the lattice laser beam 
for phase stabilization, and another one for shifting the 
phase over the range 107r with stabilization kept active. 
The short-term stability of ip is estimated to be ±0.00 77r. 
The last few optics in front of the chamber are outside 
of the active stabilization, which causes slow drift of ip 
due to the change of environment such as temperature. 
The typical phase drift is 0.057T per hour, and all mea¬ 
surements of sequential data set are finished within 20 
minutes from the last phase calibration. 


At the proper phase parameters <p x = <pz = 0 and 
ip = 7t/2, the potential depth at the center of each site 
becomes equal when Idong = Vshort = Vdiag- In this con¬ 
dition, however, the energy offset Ea becomes lower than 
Eb and Ec because of the difference in the zero-point en¬ 
ergies. We search optimal Vdi ag by single-particle band 
calculation (see also Supplementary Information for the 
derivation of Hubbard parameters). 


Band occupation measurement and sublattice 
mapping 

To measure the quasimomentum distribution of atoms, 
we turn off all the lattice potentials with a exponential 
form 

= | ^(°) exp(—4t/r) 0 < t < r ^ 

with a time constant r = 0.6 ms. The dipole trap is kept 
constant during band mapping for preventing movement 
of the trap center due to gravity and suddenly turned 
off at t = t. Because of the relatively heavy mass of 
Yb, the existence of harmonic confinement imposes se¬ 
vere restriction on the choice of mapping time r. We 
find r > 1 ms causes considerable deformation of the 
distribution whereas r > 1.5 ms is desirable to suppress 
interband transition. Due to this non-adiabaticity, up to 
20% of atoms occupying a certain Brillouin zone are de¬ 
tected in its neighboring zones, depending on the shape 
of the observed quasimomentum distribution. 
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which is utilized for the calibration of ip. 


CALIBRATION OF THE RELATIVE PHASE 

Upper images in Fig. SQJi shows time-of-flight absorp¬ 
tion images of a 174 Yb BEC loaded into an optical lattice 
at the phase ip = tv/2 and ip = 0, leading to lattice po¬ 
tential landscapes shown in the lower images. In the case 
of ip = 7t/ 2, the ground state wave function distributes 
over all three sites in a unit cell. On the other hand, 
for ip = 0 atoms are strongly localized on a square lat¬ 
tice composed of A-sites. This feature is reflected in a 
momentum distribution of a BEC, where the fraction of 
zero momentum coherent peak becomes much larger for 
ip = 7r/2 due to the delocalization. 

The phase dependence of the each momentum peak of 
a BEC is shown in Fig. £jlb- The experimental data are 
well reproduced by the non-interacting band calculation, 


TIGHT-BINDING MODEL FOR THE OPTICAL 
LIEB LATTICE 

Here we present the derivation of the tight-binding 
model describing our optical Lieb lattice. The kinetic 
energy part of the tight-binding Hamiltonian considered 
here is written in the form 

H TB = H Ueh + H BC + H s . (4) 

Here -ffueb represents the part involving nearest-neighbor 
hopping, Hbc next-nearest-neighbor hopping between 
the B and C sites, and Hs the Hamiltonian within each 
sublattice. These terms are explicitly expressed in the 
second quantized form as 




+ H.c., 


k,l 

Hbc = -JbcYI [(“ko.B^fc+bb-C + «(fe,o,s a (M-i),c) 

k,l 


H.c., 


(5) 

( 6 ) 




k,l 


Jss (“(fc,;),s“(fc-t-W),s + H.c.^ ^2 Jss + H.c.^ 


S=A,B 


+ ^2 E sa\ k ,i),s®(k,i),s 

S—A,B,C 


S=A,C 


(7) 


where is the annihilation operator on a site S 

(= A, B , C ) in a unit cell labeled by its coordinates 
(x, z) = ( kd,ld ). Inclusion of beyond-nearest-neighbor 
hopping Jbc and Jss is necessary to reproduce the band 
dispersions obtained by the first principle band calcula¬ 
tions, especially for shallow lattices. Figure 32k- shows 
a sketch of each hopping term. Regarding JbCi we con¬ 
sider terms along the direction x + z only, because our 
diagonal lattice suppresses the hopping along the x — z 
direction (see Fig. £Q]r). Similarly, hopping Jbb ( Jcc ) is 
restricted to the x- ( z -) direction, respectively. In the mo¬ 


mentum space representation a,,s = -^= e lk ' Xi - s a fc,s, 

the Hamiltonian is diagonalized with respect to the mo¬ 
mentum indices as 



T = Thieb + Tbc + 7s- (9) 

Here T is the 3x3 matrix which couples each sublattice, 
given by 


0 

TLieb = ( ~2Jcos(k x d/2) 
—2 J cos(k z d/2) 

0 0 

Tbc = | 0 0 

0 — 2Jbc cos(k x d/2 


—2Jcos(k x d/2) —2Jcos(k z d/2) 

0 0 

0 0 


0 

—2 Jbc cos(k x d/2 
k z d/2) 0 


k z d/2) 
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FIG. S 1. Phase dependence of a time-of-flight signal, a, Absorption images of a 174 Yb BEC released from the optical 
lattice with ip = 7t/ 2 (left) and ip = 0 (right), taken after 14 ms of ballistic expansion. Corresponding lattice potentials are also 
shown, b, Intensities of the coherent momentum peaks of a BEC as a function of the phase ip, measured by the time-of-flight 
experiment. Each peak is categorized into four, as shown in the right of the data plot. Solid lines show the result of single 
particle band calculation, for which horizontal offset and scale are adjusted to match the experimental data. Error bars denote 
standard deviation of three independent measurements. 


T S 


Ea — 2Jaa [cos {k x d) + cos (k z d)] 

0 Eb 

0 


0 0 

2 Jbb cos(k x d) 0 

0 E c - 2J C c cos (k z d) 


( 10 ) 


Leaving only 7Liebj we have analytic expressions for the 
eigenvalues of the Lieb lattce with hopping to the nearest 
neighbors only, as 

E± = ±2J\/cos 2 (k x d/2) + cos 2 (k z d/2), (11) 
E 0 = 0, (12) 

and the corresponding eigenfunctions 
|fc, 1st) = —= (| k,A)+ sin 9k\k, B) + cos 9k\k,C)), (13) 

V 2 

| k, 2nd) = cos0fc|fc, B) — sin0fc|fc, C), (14) 

|fc, 3rd) = —= (\k,A) — sin0k\k, B) — cos6k\k,C)), (15) 
v 2 

with tan#*; = cos(k x d/2)/ cos(k z d/2). To obtain exper¬ 
imentally relevant tight-binding parameters, we perform 
least-square fitting of the band dispersion obtained by di¬ 
agonalizing T to the lowest three bands of the first prin¬ 
ciple calculations. This procedure gives the optimal hop¬ 
ping amplitudes which reproduce the actual band struc¬ 
tures, as shown in Figures fJSb and c. Dominant contri¬ 
bution other than the nearest neighbor hopping J comes 
from Jbc and Jaa■ While the latter does not affect the 
flat band, Jbc eliminates the dark states and causes fi¬ 
nite dispersion of the flat band. The ratio Jbc/J can 
be made smaller by increasing lattice depth, or adding 
another diagonal lattice along the x + z direction. 


To derive parameters such as on-site interactions, the 
above procedure is not sufficient and we need to construct 
Wannier functions. This involves ambiguity in the def¬ 
inition of the phase of each Bloch state, and we should 
choose the phase such that resulting Wannier functions 
are well localized and minimize non-Hubbard type in¬ 
teractions. Let (u^\,u^ B ,u^ c ) (n = 1,2,3) to be the 
eigenvectors of T with the n -th energy band. The Bloch 
states obtained from the band calculations \k,n) will 
be written as superposition of the sublattice momentum 
eigenstates \k, S) = a\ s |0), in the form 

e ie ( fc ’")|fc,n)= Y, w £il*bS>. ( 16 ) 

S=A,B,C 


Here 9{k,n) is unknown phase mentioned above. We 
follow the procedure similar to that described in [l] and 
choose the phase 9 so that all Bloch states constructively 
interfere at the specific lattice site to give the localized 
Wannier function on that site: 


9(k, n) = arg 


3 -i k-Xi, s ,„,{n) ,(n 


(17) 


where ip^(x) is the Bloch wave function (x\k,n). Ac¬ 
cording to Bloch’s theorem, the above expression actu¬ 
ally does not depend on the choice of unit cell i = (fc, l). 
For each k and n, we choose appropriate reference sublat¬ 
tice S le { where the Bloch state has the largest amplitude. 
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After all, the Wanner functions Wi t s(x) can be obtained 
from the Bloch states as 

= E e " ifc ' Xi ' S E “Se" i0(fc ’ n) 4 n) (*)- (18) 

k n= 1 

Figure 33b shows the calculated Wannier functions of 
the optical Lieb lattice derived by the above procedure. 
It can be seen that the three Wannier functions wq,a> 
wo : b and u>o,c are well localized at the site A, B and 
C of the unit cell i = (0,0), respectively. It should be 
noted that, due to the difference in the confinement of 
potential wells, density at the (RC)-sublattice becomes 
higher than that of the A-sublattice. On-site interac¬ 
tions in the Hubbard model is proportional to the two- 
dimensional average density n = f dxdz|wij g(a:)| 4 of the 
Wannier functions. As shown in Fig. 33b, the sublattice 
dependence of n amounts to ~ 20 % for our lattice. 

MOMENTUM DISTRIBUTIONS IN COHERENT 
BAND TRANSFER 

In the main paper, coherent transfer of a ground state 
BEC into the 2nd band of the Lieb lattice is analyzed 
with the quasimomentum distributions. Here, we show 
a bare momentum distribution in the transfer process, 
measured by direct time-of-flight experiment. As men¬ 
tioned in the main paper, quasimomentum analysis can¬ 
not distinguish a condensate in the center of the 2 nd 
band from that in the 3rd band. In contrast, these two 
states have quite different distributions in the momen¬ 
tum space, as shown in Fig. 31b- The 2nd band is char¬ 
acterized by the prominent peaks at k = (±2 &l, 0 ) and 
(0, ±2&l), and vanishing central peak due to the destruc¬ 
tive interference of the wave function on the B- and C- 
sublattices. The 3rd band is distinguished from the other 
bands by its large population at k = (± 2 fcL,± 2 fcL)- 

Figure 31b shows the measured evolution of the mo¬ 
mentum distribution of a BEC during the coherent trans¬ 
fer scheme. The initial state dramatically changes its 
distributions and the feature of the 2 nd band is well re¬ 
produced around ~ 40 /is. At each evolution time, the 
wave function of a BEC is well described by a certain 
superposition of the three eigenstates as 

IV’) = ci 11st) + c 2 e‘“|2nd) + c 3 e I/ 3 |3rd), (19) 

where the parameters ci, c 2 , C 3 , a and /3 are chosen to be 
real and satisfy XEi c ? = 1- Since the basis states 11st), 

12nd) and |3rd) are far from the eigenstates during the 
transfer process, these parameters show complicated time 
evolutions. By regarding them as free parameters, we fit 
the experimental data with the distribution \(k\il>(t))\ 2 
to obtain the band populations c 2 (i = 1,2,3). In fit¬ 
ting, relative intensities of the coherent peaks in the re¬ 
gion —4 &l < k x < 4&l and — 4 &l < k z < 4&l are used. 


Obtained evolutions of the band populations are plot¬ 
ted in Fig. 31b, which confirms the efficient transfer into 
the 2nd band and show agreement with band mapping 
measurement in the main paper. Rather scattered data 
points might be due to the fitting to small coherent peaks 
whose intensities are largely affected by the background 
noise, or the systematic contribution from the higher en¬ 
ergy bands. 

EFFECT OF INTERACTIONS ON 
INTER-SUBLATTICE OSCILLATIONS OF A BEC 

In Fig. 4c of the main paper, we compare the mea¬ 
sured frequencies of inter-sublattice oscillations of a BEC 
with the relevant band gaps obtained by single particle 
band calculations. Systematic deviation from the pre¬ 
diction is found especially in the regime of small Sdiag- 
We experimentally examine the effect of interactions on 
the oscillation frequencies by changing total atom num¬ 
ber of a BEC, at several diagonal lattice depths shown 
in Fig. 35b- 

Assuming that the system is locally uniform, intro¬ 
duction of a local chemical potential fi(r) = n — V (r) 
leads to the expression of the total atom number N = 
/d 3 rn(/Lt(r)), where V(r) is an external harmonic con¬ 
finement and n(n) is the density averaged over a unit cell. 
For a uniform, weakly interacting BEC, the chemical po¬ 
tential has a linear dependence on atomic density, leading 
to N 2 / 5 dependence of the central density n(0). There¬ 
fore, we plot the observed oscillation frequency as a func¬ 
tion of iV 2 / 5 in Fig. 331 For each case, an extrapolation 
to the zero density seems to reproduce the expected band 
gap. Remaining errors of ~ 200 Hz can be accounted by 
uncertainty in the calibration of the lattice depths. Sys¬ 
tematic errors in the oscillation measurements can also 
come from finite lifetime (~ 1 ms) of a condensate in 
the excited bands, which matters especially in the low 
frequency modes. 

Whereas the negative shifts of the oscillation frequen¬ 
cies are observed in small Sdiag, the shifts turn into pos¬ 
itive around the Lieb lattice condition s^ag = 9.5. This 
tendency can be qualitatively explained by considering 
the density distribution of each band. For small Sdiag, Ea 
becomes much lower than Eb and Ec and the wave func¬ 
tion of the lowest band concentrates on the A-sublattice, 
whereas the 2nd and 3rd bands have population in both 
the B- and C-sublattices. The existence of a repulsive in¬ 
teraction shifts the energy upwards by a greater amount 
for the lowest band, leading to the reduction of the gaps 
to the higher bands. On the other hand, at Sdiag ~ 9.5 
the ground state wave function spreads over all sublat¬ 
tices and the amount of energy shift will decrease. The 
oscillations frequencies may depend on the initial prepa¬ 
ration, i.e., the fraction of each eigenstate. The initial 
conditions for the measurement in Fig. 351 are listed in 
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FIG. S 2. Tunneling parameters in the optical 

Lieb lattice, a, Tight binding model for the Lieb lat¬ 
tice with beyond-nearest-neighbor hopping, b, (left) Tun¬ 
neling amplitudes obtained from the tight-binding approxi¬ 
mation are shown as a function of lattice depth. The solid 
lines are linear interpolations for the successive data points. 
Signs of Jaa, Jbb ( Jcc ) are inverted to display in a log 
scale, (right) Selected diagonal lattice depth Sdiag to sat¬ 
isfy Ea = Eb = Ec . The dashed line represents Sdiag = 
siong(= s s hort) for reference, c, Band structure of the Lieb 
lattice with (siong, Sshort, Sdia g ) = (20,20,23). Dispersion re¬ 
lations of the lowest three bands of the first-principle band 
calculation (green solid) are well reproduced by the tight 
binding approximation presented here. The best fit hopping 
parameters are J = 0.282, Jbc = 0.022, Jaa = —0.026, 
Jbb = Jcc = 0.0057, Ea = 0.140 and Eb = Ec = —0.058, 
in unit of E^ 064 ^. The origin of the energy is set to the 2nd 
band at the F point. 
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FIG. S 3. Wannier functions of the optical Lieb lattice. 

a, Wannier functions of the Lieb lattice composed of the low¬ 
est three bands, at the lattice depths of (siong, s s hort, Sdiag) = 
(20, 20, 23). Scaled densities \wo,s{x, z)| 2 xd 2 are shown. Den¬ 
sity profiles of the Wannier functions on A and B sites along 
the z — 0 line are also shown in the right hand side, b, Aver¬ 
age densities of the 2D Wannier functions of the Lieb lattice 
are plotted as a function of lattice depth. The solid lines are 
linear interpolations for the successive data points. 


Smeas [siong ? S s h or t ? Sdiag] 

lst-2nd 

Smit [((S^g, Sl 0 ng), ( S Lort> S shLt)> S diag] 

|(1, meas|init)| 2 

|(2, meas|init) | 2 

|(3, meas|init)| 2 

(8, 8, 2) 

((7.6, 0), (8.7, 29), 6) 

0.964 

0.033 

< 10“ 6 

(8, 8, 9.5) 

((7.5, 0), (5.9, 29), 8.4) 

0.723 

0.275 

< 10" 6 

(8, 8, 15) 

((7.8, 0.03), (3.81, 29), 10) 

0.571 

0.426 

< 10“ 5 

1st-3rd 





(8, 8, 2) 

((8, 8), (8, 8), 20) 

0.270 

0 

0.686 

(8, 8, 9.5) 

((8, 8), (8, 8), 20) 

0.739 

0 

0.253 

(8, 8, 15) 

((8, 8), (8, 8), 20) 

0.977 

0 

0.022 


TABLE I. Initial conditions for the inter-sublattice oscillations. For each lattice depths s m eas at which oscillations are observed, 
we prepare appropriate initial state by loading a BEC into the lattice with depths Smit, and drive oscillations by suddenly 
change the depth to s mea s. Calculated overlap between the initial wave function |init) and each zero-momentum eigenstate of 
the measurement stage |n, meas) (n = 1,2,3) is also shown. 
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FIG. S 4. Momentum space observation of coher¬ 
ent band transfer, a, Numerically calculated momentum 
distributions of the zero-momentum eigenstates of the lowest 
three bands. Potential depths are set to (si ong , Sshort, Sdi ag ) = 
(8,8, 9.5). b, Evolution of the momentum distribution of a 
BEC during coherent band transfer. Each absorption image 
is taken after 14 ms of ballistic expansion, c, Oscillations of 
the band populations during coherent band transfer. Each 
population is obtained by fitting time-of-flight images with 
the momentum distribution of the superposition of the lowest 
three bands. Solid lines are theoretical expectation, using the 
parameters obtained by fitting to the band mapping measure¬ 
ment in Fig. 2c of the main paper. Error bars denote standard 
errors in fitting procedure. 
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FIG. S 5. Density dependence of oscillation frequency, a, Numerically calculated band structures of the optical Lieb 
lattice, as a function of Sdiag- Green circles are the points where the measurements in b-g are performed, with corresponding 
alphabet symbols, b-g, Density dependence of inter-sublattice oscillation frequency. Horizontal dashed lines indicate the band 
gap from single particle calculations shown in a. Error bars denote standard errors in fitting procedure. 













